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OOA would be greatly reduced. However, the Sun-S/C-Earth angle changes during the spacecraft’s orbit around the 
earth, and it changes dunnp the earth's orbit around the sun. Therefore, these sensors had to be calibrated to 
accommodate a range of Sun-S/C-Earth angles. Nonetheless, this approach reduced the error estimation effort sia- 

nifinanflv/ ® 


Celestial Attitude Reference 

The only celestial bodies which are bright enough to be seen by the MSTI 2 payload are the sun, the moon, and 
Venus. The sun could not be used because it would quickly damage the focal plane of the payload. 

The moon is bright enough to be tracked by the pay toad, but it has a significant angular extent as seen from orbit. In 
addition, its image on the MSTI 2 payload focal plane will have some unknown shape due to the cooling of the moon 
as it changes phases, and it is unclear exactly How the tracker would compute the centroid this shape. The only fea- 
sible technique to use the moon as an attitude reference would have required that the outer arc of the crescent of the 
moon be estimated, so that its center could be used as an accurate attitude reference. This implies an extensive 
development effort with high risk. 

Therefore, it was decided that Venus would be used as the celestial attitude reference. This presented its own set of 
problems. The off-axis sensitivity of the MSTI 2 payload is such that if it is pointed within 20° of the sun, any image it 
has would be washed out. Figure 2 shows the separation between Venus and the sun during the period of interest. 
Prior to Apnl 9, 1 994, Venus is within 20° of the sun. This was only a few days prior to the MSTI 2 OOA maneuvers. 


Sun-S/C-Venus Angle 



Figure 2 
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In addition this small subtended angle between the sun and Venus presented problems arising from the geometry of 
the MSTI 2 spacecraft, as will be described in the discussion on SS calibration. 


Sun Sensor Calibration 

The SS has a "square" Field Of View (FOV), nominally centered along the S/C Y axis. The sunsensorwouldrefwrt 
the position of the sun by transmitting two signed angles representing the location of the sun in the FOV . These are 
(1) SS Z, the angle between the projection or a unit vector pointing toward the sun onto the S/C X-Y plane, and the 
S/C X axis, and (2) SS X, the angle between the projection of a unit vector pointing toward the sun onto the S/C Z-Y 
plane and the S/C Z axis. 

Due to geometry constraints, it was decided that only one point in the SS FOV would be calibrated, as described 
below. A complete calibration at this single point could be characterized by a simpte bias on each of the two outputs. 
Because it is a single point, all of the non-temporal error sources could be corrected by these simple biases. 

The approach utilized to estimate these biases was purely geometric. With this approach in order to measure the 
errors with a single experiment, the full 3-axis attitude of the payload reference frame would have to be known. ^ 
However, the only hign-accuracy attitude reference available was Venus, and this was only a singl^axis attitude ref- 
erence. Therefore, multiple experiments were required to in order to make the biases observable. This purely 
geometric approach ignores the information provided by the gyros, but it also keeps the gyro alignment estimation 
and the SS calibration independent of each other. 

The approach utilized is illustrated in Figures 3 and 4. These figures are for descriptive purposes only, and 
drawn to scale. Figure 3 illustrates one procedure, and Figure 4 illustrates another procedure. TTiese two procedures 
alone are enough to algebraically determine the biases. Multiple executions of these two procedures provided statis- 
tical data and required a least-squares estimator. 

In Figure 3 (Procedure #1), the sun is located in the SS FOV at the calibration point. Also Venus is in the FOR of the 
Dayload at one location, which provides only one axis of attitude information. Using astronomical databases, Pi. the 
subtended angle between the sun and Venus, is known very accurately. Therefore it is known that the sun must lie 
somewhere along the arc shown in Figure 3. The exact location of the sun aloncj this arc is unknown because the 
Venus vector provides only one axis of attitude information. This procedure, by itself, provides only partial error 
information. In particular, it provides only e* 1f the component of the error normal to the arc in Figure 3. 

In Figure 4 (Procedure #2), the S/C is oriented so that the sun is still at the calibration point, but Venus is Jocat^d at a 
different point in the payload FOR. In effect, the S/C has slewed about the S/C-suo line frcm Procedure #1 to Proce- 
dure #2 Again this procedure, by itself, only provides partial error information, e#. Combining the results from 
Procedure #1 and #2 will provide full error information, because the errors measured in the two expenments are in 
different directions. 

The procedures depicted in Figures 3 and 4 were each executed several times in order to obtain some statistical 
characterization of the results. 

In addition to estimating the SS errors, the combination of Procedures #1 and #2 will also provide enough information 
to determine the slew angle about the S/C-Venus line for each experiment. 

In the least-squares estimator, the vector of estimated parameters was 


x = 


i 

e, 

L e 2 

Where L is the slew angle about the S/C-Venus line for procedure #i (i=1 n). e, is the bias in the SS X measurement 

at the calibration point and e z is the SS Z measurement bias at the calibration point. It was necessary to estimate the 
individual slew angles because it indirectly provided a 3-axis attitude determination for each procedure, and this is 
what is required for measuring the SS errors. 
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Figure 3 


Sun Sensor Estimation Procedure #2 



The measurement vector used for the SS estimation was: 
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m = 


ssx, 

ssz, 

ssx 2 

ssz 2 


ssx 

ssz„ 


Where SSX, and SSZ, are the two outputs of the SS for procedure #i (i=1 n). 


Because only data from the SS and the payload were used, the calibration was necessarily relative to the frame 
defined by the payload, and not relative to the S/C frame. 


Figure 5 shows the geometry of the payload FOR relative to the SS FOV. In this polar plot, the radial dimension rep- 
resents the angle from the payload Z axis, and the angular dimension represents the azimuth from the paytoadX axis. 
The solid line represents the limits of the payload FOR, and the dashed line represents the limrts of the SS FOV. Not 
all of the SS FOV is shown in Figure 5. The Reid Of Regard of the MSTI 2 payload was limited by the amount of 
travel allowed in the gimbals. 


Payload FOR and SS FOV Geometry 


so. 
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When the payload is pointed near the edge of its FOR, its image will be partially obscured by the S/C structu e, 
resulting in a dimmer image. Therefore, it is desirable to point the S/C so that Venus is as far away from the edge of 
the FOR of the payload as possible. This places a limitation on the geometnes which will simultaneously place the 

sun in the SS FOV and Venus in the payload FOR. 

Durina the time of the OOA maneuvers, the Venus was 20° to 30° away from the sun, as indicated in Rgure 2. This, 
coupled with Figure 5, indicates additional limitations on geometries which will simultaneously place Venus in the 
payload FOR and the sun in the SS FOV. 
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These limitations were the driving factors in the decision to calibrate the SS at only one point. The calibration point 
and the geometries of the SS calibration experiments, are illustrated in Figure 6. By selecting the calibration point at 
SSX = 30° and SSZ = 10°, Venus was kept as far away as possible from the edges of the payload FOR, while still 
allowing a 90° slew about the S/C-sun line between the two procedures. The small circle is tne chosen location of the 
SS calibration point, and the two Vs represent the location of Venus for the two SS Calibration maneuvers. 


Sun Sensor Maneuver Geometry 



o 


For maximum observability, the slew angle about the S/C-sun line between the two experiments should be 90°. The 
two Venus locations indicated in Figure 5 satisfy this requirement. 


Earth Sensor Calibration 


The ES is a scanning horizon sensor. It has a 2° conical field of view which sweeps out a cone with a 60° half-cone 
angle every 0.2 seconds. The axis of the ES scan cone is nominally along the S/C -Y axis. The ES measures the 
size and location of the "earth chord", which is that portion of the scan cone which intersects the earth. The outputs of 
the ES consist of <J>, the scan angle from a zero-reference scan angle to the center of the earth chord, and £1, the 
width of the earth chord. 

The basic technique to calibrate the ES is similar to the SS calibration. The primary differences are: 

(1) The geometry of a scanning horizon sensor is different from the geometry of the SS, 

(2) The ES had to be calibrated along a locus of points to accommodate a range of earth-S/C-sun angles, 

(3) No analytic horizon sensor model incorporating an oblate earth shape model was available. 

Because of the complicated geometry of a scanning horizon sensor, the location of the earth in the S/C frame was 
deduced only indirectly. Because an oblate earth shape model was required, there were no analytic expressions 
relating the location of the earth to the ES sensor readings. An iterative search was employed, using the spherical 
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earth shape approximation as a first guess. 


Because the Sun-S/C-earth angle changes during the orbit and during the year, and because only one point of the SS 
FOV was calibrated, it was necessary to calibrate a locus of ES output pairs to cover the range of possible angles. 
Each "point" in the locus corresponded to one location of the earth relative to the payload frame. 


The multiple procedure technique described above for the SS calibration was modified for the ES calibration. There 
were two sets of procedures, instead of two discrete procedures. For each procedure in the first set, the S/C was 
oriented so that the earth was at various points along the calibration locus, and all of these procedures used 
approximately the same slew angle about the S/C-earth line. For the second set of^ procedures, the earth was 
located at various points along the calibration line, but the slew angle about the S/C-earth line was approximately 90 
from the slew angle used for the first set. 


Because the ground testing included no statistical characterization, an assumed function with unknown parameters 
was used to model the ES errors. The d> and n errors were each assumed to vary as a function of the distance along 
the locus of points. Because an assumed form of the error model was used, it was unnecessary to calibrate indi- 
vidual points with pairs of parameters, as was done with the SS calibration. Instead, the error function parameteis 
were estimated collectively. It was only necessary that each individual procedure placed the earth at some location 
along the locus. It was not necessary for each procedure in the second set of procedures to place the earth at the 
same location along the locus as an procedure in the first set. 


In the least-squares estimator, the vector of estimated parameters was: 


a, 

Oa 


x = 


(Vi 

On 

C«o 

C®i 

C*z 
Coo 
Cqi 
Cq2 - 


Where ol is the slew angle about the S/C-Venus line for experiment #i (i-t n), and Cqo. c™, c^, c«o, c«,, and c M are 

the unknown parameters of the ES error model. It was necessary to estimate the individual slew angles because it 
indirectly provided a 3-axis attitude determination for each experiment, and this is what is required for measunng the 
ES errors. 


The measurement vector used was: 


m =1 


d>i 

n, 

o 2 



Where Oj and 0 ( are the two outputs of the ES for experiment #i (i=1 ,...,n). 

The least-squares estimator requires derivatives of the measurement vector elements with respect to the state vector 
elements. Because no analytic function relating O and Cl to the estimated parameters was available, there was no 
analytic expression for the related derivatives. However, the spherical earth shape model is very nearly equal to the 
oblate earth shape model, so the derivatives based on the spherical earth shape model were used, while the mea- 
surement model used the oblate earth shape model. 


Because only data from the ES and the payload were used, the calibration was necessarily relative to the frame 
defined by the payload, and not relative to the S/C frame. 


Gvro Misalignment Estimation 


310 


The ground testing of the MSTI 2 gyros was very extensive, and provided very good characterization of the ayros. 

I he gyro models, and their parameters, were determined with very good accuracy, and it was felt that no further of the 

EmSTI I S/C nee * d ^ ° nly Sit "’* ican ' em>r S0U,M ' he 9vros was ,helr 0“™** 

The sensitive axes of the one of the X-Y gyro are nominally aligned with the payload frame X- and Y-axes and the 
sensrtn/eaxes of the X-Z gyro are nominally aligned with the payload frame X- and Z-axes. The X-Lte Molma- 
tion, tor poth attitude control and telemetry, comes from only one gyro at any given time, as selected by around 
commands. The nominal mission called for using both channels of the X-Y gyro, and only the Z-axis information of 
tne x-L gyro. Therefore, it was necessary to estimate all three misalignment angles for the X-Y gyro, and only two 
misalignment angles for the X-Z gyro. y 

In order to make the gyro misalignments observable, three simple Euler-axis slews were executed by the S/C. While 
performing these slews, the payload was tracking Venus in order to make the estimated parameters observable. 
Ideally , for maximum observability, the three slew axes should be mutually orthogonal. However, the limited extent of 
the payload FOR constrained the slew axes to be non-orthogonal. The geometry of the three slews is illustrated in 

°£ en,at . 1 ?? ?, f the . a []9ular velocity vector for each slew is indicated by the small "o"s, and the dashed 
circles indicate the path followed by Venus through the payload FOR. 


Geometry of Gyro Misalignment Estimation Maneuvers 



estimation, it was assumed that the two axes of each gyro are perfectly perpendicular, 
that the misalignment of each gyro was independent of the misalignment of the other gyro. 


It was also assumed 


'i!^ was comman ded to perform the three slew maneuvers, it complied by using the misaligned gyros. 
Therefore, it was necessary to also estimate the parameters of the Euler-axis slew as well as the gyro misalignments 
I he estimated parameters were: a 


(1) The orientation and magnitude of each slew axis in the payload frame coordinates, 

(2) The orientation of each slew axis in inertial space, 


311 


(3) The initial slew angle of each slew maneuver, 

( 4 ) The X, Y, and Z angular misalignments of the X-Y gyro relative to the payload frame, and 


( 5 ) The X and Y angular misalignments of the X-Z gyro relative to the payload frame. 


it wac noressarv to exDress the equations describing the slew maneuvers in terms of the orientation of the slew axes 
h taSSrnEStf (X to (orcelffiSuniRfl slew girametera to apply to an Euler-axis slew toed in mental space. It 
was necessaty lo estimate these orientations because they were not directly measurable. 


In the least-squares estimator, the vector of estimated parameters was: 
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x = I DecL. 
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Where 6 a, is the initial slew angle ol the SIC at the beginning erf ^ nranawer #i 0=1 °" ** 

maneuver #i (i=1 ,2,3), are the three mislaignment angles of the X-Y gyro (r-1 ,2,3), and y* are two or me mree m.s 

alignment angles of the the X-Z gyro (i-1 ,2). 

The X axis data from the X-Z gyro was not available, so the third misalignment angle of the X-Z gyro, Y 23 . ' s neither 
required nor observable. 


The measurement vector used was: 
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m = 


v P „ 

V p21 

Vp,i 

A 

tain 

A 

a* 2i 

A 

ffci 

Vp,2 

Vp22 

V p» 

A 

t*W 

A 

^22 

A 

t^32 


v pln 

^p2n 

Vpv, 

A 

ta,i„ 

A 

tai2n 

I A 
L «*3n 


VVhere v* is the component of the Venus unit vector along the p, direction at time t (i=1 ,2,3; j=1 n), cL is the angu- 

lar velocrty measured by the X axis of the X-y gyro at time tj (i=1 ,—,n), co^ is the angular velocity measured by the Y 
axis of the X-Y gyro at time tj (i=1 ,...,n), and ofe* is the angular velocity measured by the Z axis of the X-Z gyro at time tj 

\* * ,...|lly. 


Conclusion 


Most of the enor sources in the attitude knowledge of the MSTI 2 spacecraft were removed analytically by using the 
payload as a high fidelity attitude reference, and using this information to deduce the enors in the ACS sensors. This 
determined the errors relative to the payload, and eliminated the need to align the sensors relative to the spacecraft. 
Those errors which could not be eliminated analytically were avoided by careful mission design. 
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